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I Abstract 

In this paper, a link between g-difference equations, Jacobi operators and orthogonal poly- 
nomials is given. Replacing the variable x by g"" in a Sturm-Liouville g-difference equation 
we discovered the Jacobi operator. With appropriate initial conditions, the eigenfunctions 
of such operators are either q-orthogonal polynomials or the modified g-Bessel function and 
a newborn the g-Macdonald ones. The new Polynomial sequence we found is related to the 
g-Lommel polynomials introduced by Koelink and other [9, 10, 11]. Adapting E. C. Titch- 
. marsh's theory, we showed the existence of a solution square-integrable only in the complex 

case. As application in the real case we gave the behavior at infinity for q-Macdonald's 
function. Finally, we pointed out that the method described in our paper can be generalized 
1^ I to study the orthogonal polynomial sequence introduced by Al-Salam and Ismail [2]. 
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^ '■ 1 Introduction 

p_ 

In this paper we study the relationship between three principal topics, namely 
1. g-difference equations 
^ ■ 2. Jacobi operators 

3. Orthogonal polynomials 

Each of these subjects has been intensively studied by various authors and some links have been 
^ I established to joint them. For example, the theory of g-orthogonal polynomials was built while 

■ studying (^-difference equations. The spectral study of Jacobi operators was started by many 

authors but here we showed their relationship with the g-difFerence equations. So what is new 
in this work? 

There are two ways of seeing a g-difference equation: 

- The first is to study a g-equation with one variable x of the Sturm-Liouville kind and find 
power series solutions. 

- The second rises from a simple observation, which turns out successful, that is to replace in 
the g-difference equation the variable x by which brings up the Jacobi operator. With 
appropriate initial conditions, the eigenfunctions of such operators are either g-orthogonal poly- 
nomials or special functions. This, immediately, able us to highlight the link between g-special 
functions and g-orthogonal polynomials. One of the consequences of this connection is either 
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their generating function or their measure of orthogonahty. In our case, the special functions 
are the modified g-Bessel and among them we have a newborn the g-Macdonald function which 
we investigated its properties more closely with the great help of the g-Bessel Fourier analysis 
tools [3]. 

It turns out that this g-Macdonald function is different from the g-Macdonald function of third 

kind introduced by Rogov [12]. On the other hand we found a new family which is related to the 

g-Lommel polynomials introduced by Koelink and other [9, 10, 11]. Thus we could determine 

their generating function and their orthogonality measure (which is discrete). 

To enrich our paper we applied the methods of the two classical theories in our study. 

In the case of the Sturm-Liouville equation we were able to adapt E. C. Titchmarsh's theory 

[15, 16] to show the existence of a solution square- integrable only in the complex case. 

In the real case we found another approach giving, as application, the behavior at infinity of the 

g-Macdonald's function. 

In the case of the Jacobi operator we used the tools of moments theory to obtain some additional 
information on the subject. 

We completed this work by noting that the method described in our paper can be generalized 
to study the orthogonal polynomial sequence introduced by Al-Salam and Ismail [2]. 



2 Preliminaries and g-Bessel Fourier transform 



Throughout this paper, we consider < g < 1 and v > —1 such that ^ N. We adopt the 
standard conventional notations of [5]. We put = {g", n G Z} and for a complex number 
a 



n-l 



and 



(a; g)o = 1, (a; q)n = Y\.(^ ~ ^^^)^ ^ = +°°' 

1=0 

{q, q)n 



n 
m 



n—m 

The g-hypergeometric series (or basic hypergeometric) series r<^s is defined by 



ai , . . . , Or 
bi,...,bs 



y h i^-^^Ubi;q)k---{bs;q)k ' ^ ) 



1+s— r 



whenever the series converges. 

The third Jackson g-Bessel function Jy (also called Hahn-Exton g-Bessel functions) is defined 
by the power series [16] 



T / \ {.q ^ ') q)oo I, 

Jv\x\q) = — X 

{q;q)oo 



EM) 



n(n+l) 

g 2 



„2n 



n=0 



{q'"^^;q)niq;q)n 



and has the normalized form 

oo 

Mx;q) = ^(-1) 



n(n+l) 



„2n 



n=0 



iq,q)niq''+\q)r 



The g-derivative of a function / is defined for ,t 7^ by 

f{x) - f{qx) 



(1 - q)x 
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and the g-Bessel operator by 



1 



The g-Wronskian was introduced in [4] as foUows 

WM, h) = q-\l - qf [k-^D^f{x)h{x) - A-^D^h{x)f{x) 
where 

Ag/(x) = f{qx). 

The Jackson's g-integrals from to a and from to oo are defined by [7] 

pa °° foo °° 

/ /(xKx=(l-g)aJ]g"/(ag"), / f{x)d^x = {\-q) J] q^f{q^) 
Jo n=0 Jo n=-oo 

For 1 < p < oo, the space Cq^p^n denotes the set of real functions on M+ for which 



\f{x)\Px^^+^d,x 



1/p 



is finite. 



The g-Bessel Fourier transform Fq^u is defined by [4, 8] 

/•oo 

J'q,vf{x) = Cq,, / f{t)j,ixt, q^)t^''+^dqt, G M+ 

Jo 

1 (g''^+';«2)oo 



where 



1-q (g2;g2)^ 



In [3] the following results was proved: 
Proposition 1 For all v G R/Z_ we have 

W,q^)\< 



n 2,, ^ {-e-.q^Ui-q^'^^-A r 1 



if n > 

^n'-{2u+l)n if n<0 



Proposition 2 Let n,m E 7. and n ^ m, then we have 

/•oo -2n{u+l) 
4,. / Mq'^x, q^)Mq"'x, q^)x'-+UqX = \ 5nm. 

Jo ^ — q 



In particular 



/•oo 

/ Mq''x,q^)x'''+^dqx = 0. 
Jo 



Theorem 1 For f G jOq^i^^ we have 



For f,Fq,y{f) G Cq^i^y we have 



:FlAmx) = f{x), VXGM+ 



\\^q,v{f)\\q,2,v - \\f\\q,2,v 



We end this section by the following result (see [4]): 
Remark 1 Let f,g& ^g,2,u such that Ag^^^f, Aq^^g £ Cq^2.u- If we have 

Dgf{x) = 0(a;-^) and Dqg{x) = ©(a;"'') 

as X \.Q then 

{Aq,,f,g) = {f,Ag,,g), 
where {.,.) denotes the inner product on the Hilbert space C.q,2,v 

3 Modified g-Bessel functions 

Definition 1 We introduce the modified q-Bessel functions as follows 

Iy{x, q^) = jy{ix, q^), = -1, 



i^y{x,q^)=x '^"j-uiq "ix^q 



and 



K,{x, q^) = c,,, / [1 + t"] Utx, q^)t''+^dqt. 
Jo 

This function can be viewed as a q-version of the Macdonald function (see [14, p.434j)- 
Theorem 2 The q-Macdonald function K^, G Cg^i^i, and we have 

J^,,,(K,)(x) = [1 + x^] ~\ Va; G M+ 
Proof. To show K,y G jCq^i^i/ we consider three cases: 
i. > 0; ii. — 1 < < and ill. = 0. 
Case i. For x < 1 we have 



(1) 



poo 

x'''+^\K,{x,q')\ < Cg,J 

Jo 

roo 

Jo 



21 



X 



\Ut,q')\t^''-'dgt}x' 



< <c, 



f)\t'^-'dgt[x' 
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and then / \Ki,{x,q'^)\x'^'''^^dqX < oo. On the other hand, for a; > 1 we have 
Jo 



Jo 



= -Cg,, [ ^ Ag,,j,{xt,q^)f''+^dgt 

Jo 

= -Cq^u Ag,^[l + t2] j^{xt,q^)t^''+^dqt (from Remark 1) 
Jo 

roc 

= -Cq,u / Uit)j,{xt,q^)t'"'+^dqt 

Jo 



where u{t) = Ag^,^ [l + 1^] ""^ is a bounded function on R+. Multiplying x'^K,^{x,q'^) by 
and making an upper bound for the quantity we obtain 



x'''+'\K4x,q' 



poo 

< Cq,, / \u{t/x)\\Mt,q^)\t'"'+^dqt 

Jo 

< ^Cq,4u\\q,oo |> K'"+'c/,t| , 

/oo 
\Ki,{x,q^)\x^'^'^^dqX < oo. 

Case ii. In fact lim Ki,{x, q^) exists and then / \Ki,{x)\x'^^^^ dqX < oo. Using the same method 

/oo 
\Kiy(x, q'^)\x'^'^^^ dqX < oo. 



Case iii. For x < 1 we have 



x'^\Ko{x,q'^)\ < Cqfi 

Jo 



21-1 



1 + ^ 



x^ 



\joit,q^)\tdqt 



< < C. 



■5,0 



X 



1 + - 



\jo{t,q^)\dqt } X 



< |cg,o J \joit,q'^)\dqt^ X 
and then / \Ko(x,q'^)\x'^^~^^dqX < oo. The proof is identical to case ii. 

^0 

This proves that G J0.q^i^,y and then (1) is a simple consequence of the inversion formula in 
Theorem 1 ■ 

Proposition 3 Let f and h two linearly independent solutions of the following q- difference 
equation 

Aq^yy{x) = ± y{x). 
Then there exists a constant c{f, h) ^ such that 



X 



2u 



f{x)h{qx)-f{qx)h{x) =c{f,g), Vx G M+ 



Proof. With 



0, 



Dq [y ^ W^(/, h)\ (x) = [Ag,J{x)h{x) - h{x)Ag,,fix)\x'"'+' = 

there exists a constant c such that 

- qfx^''+^ [A-^DJ{x)h{x) - A-^Dgh{x)f{x)\ = c, G M+ 

To completes the proof, it suffice to use the definition of the operator Dq. m 

Proposition 4 The functions x i-)- ju{^x,q^) and x i-^ 'y,y{Xx,q'^) are two linearly independent 
solutions of the following equation 



\,vy{x) = -X^y{x). 



In addition 



c{j.,l.) = {q-^'-l). 
Proof. For all a G M such that a 7^ —1, —2, ... we have [3] 

Aq,aJa(Ax,g^) = - \^ j a{\x , q^) . 

Let ji = q~^X, then we have 

-2v g""i-.(g-^-"Ax, q^) - (1 + q^-)j^,{q--Xx, q^) + j-.(g^-"Aa;, q^) 

x^ 

-2.i-.(r'-"Aa;, g2) _ (1 + g-2-)j_,(g--Ax, q^) + q-''^ j-M'^'^x, q^) 



\,vlv{Xx,q^ 



X 



2v - 

q X 



x^ 



g2u^-2u j-AQ ^t^x, q^) - (1 + g ^•')j-^{iix,q'^) + q '^''j-^{qnx,q'^) 



x^ 



2u-2u 



q"x " [-iJ?j-^{jix,q^ 



= -A^7y(Ax,g^). 
On the other hand we have 

c{ju,lv) = q~'^'"ju{x,q'^)j-u{q^~''x,q'^) - jt,{qx,q'^)j-t,{q~''x,q'^). 
When .T — 7- we obtain the result. ■ 

Proposition 5 The functions x 1— I^{Xx,q^) and x Ki,{Xx,q^) are two linearly independent 
solutions of the following equation 



\,vy{x) = X^y{x). 
Proof. Let gx{x) = X^^''+^^ K^{Xx,q'^) then 



(2) 



gx{x) = Cq^y I 

Jo 



1 + 



A2 



Utx,q^)t"'+Hqt. 



A2 



9\{x) 



f 

Jo 



2 1-1 



A2 



Mtx,q^)t^''+%t 



roo 

= Cq^u / ju{tX,q^)t^^+^dqt = Q. 

Jo 



Note that 



1 - 



A2 



1 + 



A2 



The function G C'q,i,v but 7,^ ^ ^q,i,i^- Hence, we conclude that they provide two hnearly 
independent solutions. ■ 

Lemma 1 Let A G C such that A ^ U g^'M^ i/ten w;e have 



>(g-"A,g2) (g-2.+2^g2)^ (^-2^g2)^(^2^2^^2)^ • 

Proof. Let x G C*\M+ then we have the following asymptotic expansion as \x\ oo 



(3) 



jaix,q ) 



2x {x^q^,q^' 



If n — >■ oo we obtain 



■ / -n\ 2\ W ,q )oo 

Ju{q A,g ) - 



(9 



'2^+2, g2)<^ 



VA ^ M+ 



On the other hand we have 

(g^-^"A2,g2)oo = (-l)"A%-(-i)(A-2,g2)^ (g2^2^^2) 



oo 



Hence, when n — >■ oo 

. , -n, 2, (-l)"A2»g-"("-l)(A-2,g2)^(g2^2^g2) 

jj.(g A,gj~ 

and 

J 

This implies 



VA^ 



. -n-., „2^ (-l)"g-^-"A2''g-»(»-l)(g2^A-2,g2)^(g2-2.^2^g2^^ 

,/(g A,g j (g-2i.+2^g2^^ ' vA^gM^. 



2.n J-^(g"""'^A,g2] 

g ~ 



g^)oo (g^'^A-2, g2)„(g2-2^A2, g2)^ 



>(g-A,g2) ~ (g-2-+2,g2)^ (^"2, g2)„(Q2A2, ^2)^ 



, VA ^ R+ U g' 



Hence, 

2,„ j-.(g— A,g2) ^ (g^-+^g^)oo (g2-A-2, g2)^(g2-2.A2^ ^2^^ 

nToo'^ j.{q---\q^) (g-2-+2,g2)^ (A-2^g2)^(g2p^g2)^ 

This completes the proof. ■ 
Proposition 6 We have 



-, VA^M+Ugl 



where 



and 



a,, = < 



Kt,{x, q^) = a„ 7T^{x, q^) - Pulu{x, q^) 

(g^^+^g^)oo (-g^^g^)oo(-g^~^^g^)oo 
(g-2'^+2,g2)oo (-I,g2)oo(-g2,g2)oo ' 

f if i^>0 

/•oo 

/ [l + t'l-^'^^+^qt if u<Q 



(4) 
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Proof. The functions x i— )■ /jy(Ax, q^) and x i— )■ 'Ky{\x, q^) arc two linearly independent solutions 
of (2). Then there exist two constants a^, and jin such that (4) hold true. Now we can write 



On the other hand 



Using Theorem 2 we see that 
Then it is necessary that 



lim '^,q'^) = +00 

n— >-oo 



lim K^(g-",g2) = 0. 
n— >oo 



lim 

n— >oo 



/9. 



= 0. 



Formula (3) with X = i lead to the result. To evaluate a^, we consider two cases: 
If 1/ > then 

roo 

a,y = lim x'^''K,y{x,q'^) = „ I ji,{t,q^)t^''~^d„t. 
Using an identity established in [8] 



/>oo 

/ t-^Je{q^t-q^)J^{q'^t-q^)d,t 
Jo 



(1 - q'jqM>^-'^)+(m-n)n i'i 



1+X+e-ti „2/i+2.„2 



l-A+At+6» 1-A+M- 



2. 2m-2n+l+A+6»-/i 

y ^ q 



valid for ra, m G Z and ^, /Lt, A G M such that 

l-A + 6' + /x>0, 2m-2n + l + A + 6'-/x>0, 
we conclude that (0 = z^, /x = 0, A = 1 — i^, ra = 0, m ^ 00) 

{q' 



,2 2x 

) q joo 



{q^-,q^ 



If -1< 1/ < then 



"i/ = --5- lim Kv{x,q^) 



f2i/+l 



dnt. 



Corollary 1 ^ds a direct consequence, we have 

c(I,,ir,) = a,(g-2--l). 
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4 Jacobi operator and ^-difference equation 



A tridiagonal matrix of the form 



J = 



ao 

















h 


ai 











ai 


62 


02 











a2 


bs 


03 






•A 



is a Jacobi operator, or an infinite Jacobi matrix, if 6, G M and a, > 0. 

We consider J as an operator defined on the Hilbert space ^^(N). So with respect to the standard 
orthonormal basis {efejfeeN of ^^(N) the Jacobi operator is defined as 



J Ck 



akek+i + bkCk + ak-iek-i, k>l, 
aoei + boco, k = 0. 



Associated with the Jacobi matrix J there is a three-term recurrence relation 

(A) = AP„(A), n > : 



(5) 



If we take the initial conditions P-i{x) = and Po{x) = 1 we obtain a sequence of orthogonal 
polynomials (Pn)- We assume that P„ is of degree n with positive leading coefficient. 

The polynomials of the second kind {Qn) are generated by the three-term recurrence relation 
(5) but with initial conditions Q-i{x) = —1 and Qo{x) = 0. Consequently, (P„) and {Qn) are 
linearly independent solutions of (5) and together they span the solution space. Notice that 
Qn{x) is a polynomial of degree n — 1. The orthonormal polynomials {Pn) and the polynomials 
of the second kind {Qn) play a crucial role for the theory of the moment problem. 

If {Pn) satisfy the above three-term recurrence relation (including the initial conditions) for some 
real sequences (a„) and (6„) with a„ > 0, then it follows by Favard's theorem that there exists 
a positive measure /it on M such that the polynomials {Pn) are orthonormal with respect to /x. 
The recurrence coefficients in (5) contain useful information about the moment problem. Car- 
leman proved in 1922 that the moment problem is determinate if 



00 ^ 



n=0 



00. 



(6) 



If ji is determinate, then Y^^=o \Pn{z)\^ = 00 for all z G C\M and all 2; G M except for at most 
countably many points z where fi{{z}) > 0, in the case of which 



\n=0 / 

If fj, is indeterminate, then 

00 00 

5^ |Pn(z)|2 < 00, <oo, VzGC. 



(7) 



n=0 



n=0 
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There exists close relationship between g-difference and Jacobi operators. Let us describe briefly 
this relationship. 

Given an eigenfunction / of the following g-difference operator 

Hg,,f{x) = A,,,/(x) - V{x)f{x) = -Xf{x), Va; G M+, 

then we have 



1 



x^ 



1 + 



■2w 



+ V{x) 



2u 



fix) + ^f{qx) = -Xf{x). 



Assume that f{x) = x ^'^~^^^g{x), then 



x^ 



-g{q ^x) 



1 + q"^" 
x^ 



q 

g{x) + -^g{qx) = -\g{x). 

X'^ 



(8) 



If we replace xhy q and g{q ^) hy gn in (8) then we obtain 

Jgn = CLngn+l + ^n^n + On-lgn-l = — Ag'n, 

where 

an = bn = + q'n - V{q-^). 



(9) 



Note that: 

• {dn} is an eigenfunction of the Jacobi operator J if and only if / is an eigenfunction of the 
^'-difference operator Hg ,^. 

poo 

■ {gn} G ^^(N) if and only if / |/(x)px^''+^dgX < oo for some constant a G M+. 

J a 
poo 

. {gn} e e^{Z) if and only if / \f{x)fx'^''^^dgX < oo. 

Jo 

Theorem 3 We assume that X is a complex number other than real values. There exists a 
unique solution f (up to constant factor) of the q-difference equation 



satisfying 



= [V{x) - X]f{x), 

POO 

•J a 



(10) 



x)\^x^''+^d„x < oo, 



for some constant a G . 

Proof. If F{x) = F{x, A) and G{x) = G{x, A') satisfy (10) then, 



(A' -A) / F{x)G{x)x'^''+^dgX 

J a 

F{x)[v{x)G{x) - Ag,,G(x)} - G{x)[v{x)F{x) - A,,,F(x)}] x^^+V.x 
= -J^ [f(x)A,,,G(x) - G{x)Ag,,F{x)\x'"'+%x = a''^+^WaiF,G) - b''^+^Wb{F,G). 
where Wa{F,G) is the g-wronskian. 
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If A = /X + iry, A' = A and G = F, then we get 



2u / \F{x)\'^ x'^'^+^dqx = ia'^''+^Wa{F,F) - ib^'^+^WbiF, F). 

J a 

In the sequel we assume that rj > 0. 

Now let (p{x) = 4>{x, A) and 9{x) = 6{x, A) be two solutions of (10) such that 

J 0(a) = esina, Dq(p{a) = — ecosa 

I ^(a) = ecosa, Dqd (a) = e sin a ' 

where a is real and 



It follows that 
and 

Therefore 



b'^''+^Wb{(l), 9) = a^^+^Waiy^, 9) = sin^ a + cos^ a = l, 
a2-+iH^„(0,0) = a^''+'Wa{9,9) = 0. 



a'^''+^Wa{e + l(j),9 + = I - J = 2ilml. 
Let z be a complex variable such that 



I = /(A) 



9ib)z + Ag^Dg9{b) 

' 4>{b)z + Aq^Dg(t>{b)' 



If z describes the real axis D = {z G C, Im z = 0} then I describes a circle in the complex 
plane. 



Here / = oo correspond to z 
therefore 



A-^i^,0(6) 

m ■ 
h = - 



Also 



Im 



-I 



The center of correspond to the conjugate, 

Wb{9,^) 

A-iL>,0(6) 4>-'Dg<f>{b) 



m 

b 



2 \m? b^^+'imi'Ja 



\(l){x)\^x'^''+^dqX > 0. 



The exterior of Cb correspond to Imz > 0. Now I is inside if Imz < 

i.e. if _ _ _ 

lA-^D,m + A-^Dg9{b)_ ^ lA-^DgcPib) + A-^Dg9ib) 



l<j){b) + 9{b) 



l<P{b) + 9{b) 



>0, 



equivalently 



Wbi4>, 0) + lWbi4>, 0) + lWbi9, 4>) + Wb{9, 9) 



>o, 
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that is 
equivalently 



iWb{e + i^,e + let)) > 0, 

lm.1 



[ \0 + l^\ 

J a 



Since — - 



is on Cb (for z = 0) the radius of Cb is 



n 



A-^Dqe{b) Wh{e,<p) 



2riJ^\<j){x)\^x^''+^dgX 



lib < b' then C^/ includes Cb (r^/ < r^). As 6 ^ oo, the circle Cb converge either to a hmit-circle 
or to a limit-point. 

If m = m(A) is the limit-point, or any point on the limit-circle, 



/ \e + m(f)\ 

J a 



for all values of b. Hence 



POO 

J a 



Imm 
Imm 



Note that if Im(A) = < then it suffices to consider the following ^'-difference equation 

A,,J(x)= [F(x)-a]7(x). 

Using (9) and (6) we see that the moment problem is determinate. By (7) we conclude that 
when the moment problem is indeterminate we are in the limit-circle at infinity. So in our case 
we are in the limit-point at infinity. ■ 

In the following we assume that A is a real value and there exist a G and (5 > such that 

V{x) - A > 5, G R+, x>a. (11) 

We denote by 

/n = /(g-"), yn = q-''4^, if /n-17^0. 

Jn-l 

Lemma 2 Assume that f{x) 7^ for x big enough and 



roo 
/ 1/(^)1 



2 2i/+l J ^ 
X ^ dqX < 00, 



(12) 



then when n is big enough we have 
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Proof. In fact ii a = q "° then 



n=no 



which impUes 



Hm sup 

n— >oo 



-(i.+l)(n+l)|j(-^-n-l 



< 1. 



When n is big enough we have 



-(!^+l)(n+l)|j(^^-n-l-)| 



1 

< -, 



and the result follows. ■ 

Theorem 4 Lei f be a solution of the q- difference equation 

\,uf{x) = [V{x) - A] fix), Vx G M+ 

r/ten either lim /(x) = ±00 or / satisfies the conditions of Lemma 2 and then when n big 
enough we have: 

i. f „ > and lim z/„ = 0. 

n— >oo 



ii. There exist c, a > suc/i i/iai 



l/nl < oc^q'' 



Proof. We have 



Let a G such that 



Assume that /(o) > 0. 
★ If Dqf{a) > then 

which implies 

This gives 



g(l-g)-V'^+M^(^)-A]/(x) = g(l-g)-V'^+^A,,,/(x) 

= D,[y^y'-+'DJ{q-'y) 

V{x) -X>6, Vx > a. 

DJiq-^'a) > DJ{a) > 0, 
/(g-"a) > {q--a)DJ{a) > 0. 



lim fiq-"") = +00. 
n— >+oo 



* If Dqf{a) < then there are three possibilities: 
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• The function f{x) changes sign, it can do so at most once and then goes to — oo. (replace 

/ by — / and use previous argument). 

• The function f(x) remains positive and Dqf{x) change sign then f{x) goes to +00. 

• The function f{x) remains positive and Dqf{x) remains negative for x > a, then / 

decreases steadily, and then goes to a positive limit. Hence Dqf{x) — ^ 0. Moreover 

f Aq,,fix)x'"'+%X = f [V{x)-\]f{x)x'''+^dqX 

< -q-\l-qfa^''+^DJ{q-^a), V6 > a. 

Hence 

6 / f{x)x^''+UqX < / [V{x) - X] f{x)x^''+'dqx < -q-\l - qfa^''+'Dqf{q-'a), 
J a J a 

which proves that / satisfies (12). 
Similar arguments hold if /(a) < 0. 

Assume that / satisfies the conditions of Lemma 2. We have 

f{q-'x) - (1 + q-'nfix) + q^^'fiqx) = x^ [q{x) - A] f{x). 



Let 

then we have 



9n = V{q-^), 



U+l - (1 + q'ntn + = 9-'" [qn " A] /„ 

9-'"[gn-A] + (l + g2-) 



—i/fn+1 , vfn—1 

q — ? — \-q —r- = q 



Set 



we get 



Prom Lemma 2 we see that 



fn fn 



g-'" [qn-X] + {l + qn 



^n+l = An . 



A-n 



|fn+l| < 1- 



As V{x) satisfies the condition (11), we obtain 



Hence 



lim An = +00. 

n->oo 



z/„ > 
lim = 
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On the other hand 



Hence, as n tends to +00 



< (g'^ + g-'^ + l). 



1 



l^n [qn - A] 

Then when n is big enough we have 



< -^^^ —. ^—q'^"- 0. 



2n+v 



[qn - A] ^ 2 

r r.'if 

\fn\ < 



9'"|/n-l|, 



which leads to the result. ■ 

Corollary 2 The q-Mackdonal function satisfies the following properties when n is big enough 
i. The sequence Ki,{q~'",q^) have a constant sign and 



ii. There exist c,a > such that 



K,{q~'',q^)\<ud'q'' 



Proof. From Theorem 2 we see that the g-Mackdonal function satisfies the condition of Lemma 
2. Then we apply Theorem 4 with F = and A = — 1. This give the result. ■ 

5 Orthogonal polynomial 

Let y = in (9) and consider the three-term recurrence relation 

g2»+-+ip„+i(A) - q^^il + q^-)Pn{\) + g2-+-ip„_i(A) = AP„(A). (13) 

with the initial condition P_i = and Pq = 1. Using (6) we see that the moment problem is 
determinate and so, there exists a unique positive measure /x on M such that the polynomials 
{Pn) are orthonormal with respect to fi. 

Remark 2 The q-Lommel polynomials rn{-,w,q) are given as follows [9] 

\w\~^ q~^~'^rn+iix,w,q) + + w~'^)rn{x,w,q) + \w\~^ q^~'^rn-i{x,q,w) = xrn{x,w,q) 

with initial conditions 



where 



r-i{x, w, q) = 0, ro{x, w,q) = 1 
0<q <l, > 1. 



(14) 
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So it is easy to show that if w = q ^ then we obtain 

Pn{x) = rn{x,w,q~^). 

Here Pn{x) are the polynomials satisfying (13) and are not a particular case of the q-Lommel 
polynomials because they do not satisfy the conditions (14)- We can consider Pn{x) as the 
q^^ -Lommel polynomials. 

Moreover if we replace x by q"^ and g{q^) by Qn in (8) then the recurrence relation (13) becomes 

So with the choice w = q~'" and the condition v > Q we obtain Pn{x) = r„(a;, lu, g^) which are 

a particular case of the q-Lommel polynomials. The polynomial sequence (13) is not q-classical, 
in fact, from [6, p. 65-67], we can see that the recurrence coefficients are neither classical nor 
shifted- classical. 

Proposition 7 Let f and h be two linearly independent solutions of the following q-difference 
equation 

Aq,,y{x) = Xy{x). (15) 

Then 



Pn{\) = q 



-n(i/+l) 



f{q-^)h{q) - h{q--)f{q) 
f{l)h{q)-h{l)f{q) 



Proof. The polynomials (Pn) and {Qn) are linearly independent solutions of (13) and together 
they form the space of solutions. Then there exist two functions and m such that 



We have 
This implies 

For n = — 1 we get 
and for n = 
Then 



gn = e{\)Pn{X) + m{X)Qn{\). 
f{q~n = q""^"^'^ kA)Pn(A) + m(A)Q„(A)' 



m(A) 



= /(I)- 



,n(v+l) 



f{q-n = q 

Also for the function h we have 

which give the result. ■ 

Corollary 3 We have the following identity 



/(l)P„(A)-g(''+i)/(g)g„(A) 
h{l)Pn{\)-q^''^''^h{q)Qn{\) 



(16) 



a) >(g "A,g2 



n{i/+l) 



MX, q^)Pn{-X^) - q^'+'huiqX, q^)Qn{-X^) 
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c) I.(g-A, g2) = Jj^^^^ q^)P^{X^) - g(-+i)/,(gA)g„(A2 

d) K,{q-^\ q^) = [k,(A, q^)Pn{,\^) - q^^+^^K,{qX, q^)Qn{\^) 



f) P„(A2) = g-^^'^+i) 



A2- 



A- 



2u 



a,(g-2- - 1) 



Proof. The identity a) and b) are simple consequences of (16) if we replace A by — A^ in equation 
(15). For the identity c) and d) we replace A by A^. To prove e) and f) we use Propositions 3 
and 7. ■ 



Corollary 4 We have 



^ m=0 



where 



'"m,n 



q ^ 2<Pi 



-2m -2m+2v 
r,2v+2 



_ ^2n{m+v) 



—2m —2m+2v 

q ) y 

„2u+2 



g2^^4+2(.-n) 



Proof. Using the identities e) in Corollary 3 we obtain 



Pn(-A' 



2x q 



-n(i/+l) 



-n{u+l) 



A 



2v 



Uq-^K q')lu{qK q') - >(gA, q'h^iq-'^K q') 



q-'lAq-''^ q'jj-Aq'-'^ q') - ^'"^^(^a, g^)j_.(g-"-^A, q 



n—u \ „2^ 



Taking into account formula [16, (6.4.4)] for the product of two Hahn-Exton g-Bessel functions 

2' 



°° qm(m+l)y2m 
ja{ax,q^)jp{hx,q^) = ^(-1)"" , o on / or^o ox 291 



m=0 



{q^q^Uiq^^'^^q^) 



-2m -2m-2/3 

q t q 



q 



2a+2 



2m 



we obtain 



°° m{m+l)+2m{l~u)-2u 

q-"'jAq~''x,q')j-u{q'-''x,q') = ^(-1)"^ .^. . _.„^. 2<Pi 

m=0 



{q^q^)m{q-"'+^,q^)r 



-2m -2m->r2u 

1 q 

„2i/+2 



q2^g2{u-n) 



and 



g2-j,(gx,g2)_j-_,(g«-^,g2) ^ £(-1)^ 



m=0 



gm{m-'rl)-i-2m{n—u)+2nu 



291 



-2m -2m+2i/ 

1 q 



q 



2i/+2 



^2^^4+2(.-n) 



which give the expression of the coefficients am,n- 
Note that am,n = if m > ra. ■ 
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Proposition 8 The generating function for the polynomials Pn is 
^ 1 



{1 - q''-H){l - q-i'^+^H) 



292 



n=0 

As a consequence we have 

n 
m=0 



q^;-q-(^+-)xt 



, q^-H and q-^^+^h 



n 






2(t>l 


m 









2m+2 2(m-n) 
„-2n 



q2.^2{u-m) 



Proof. As in [10, p. 13] we will derive a generating function for the polynomials P„. Multiply 
(13) with and sum from n = to oo we have 

oo oo oo oo 

n=0 n=0 n=0 

Let us introduce the generating function 



n=0 



then one gets 



With G{x, 0) = 1 we obtain 



G{x,qH) 



G{x,t) = Yt^Pn{x) 
n=0 

l + q^-''xtG{x,t) 



(1 - q^+'^t) (1 - q^-'^t) ' 



qm{m—l)—{l+i')m 



{xty 



m+l 



{1 - q''-H){l - q-(''+^H) 
Using the (^-binomial theorem one obtains 

1 



2'?!'2 



q^. -q-i^+-)xt 



1 



= l(t>0 



' q2m+2 















{q-(^+^H,q'^) 



m+l 



J2m+2 



, \q''-'t\<l. 
, \q-^^+'h\<l. 



This yields 
G{x,t) 



^m(m-l)-(l+i/)m+i(i/-l)-fc(i/+l) 

m,j,k=0 



( 2m+2 2\ / 2m+2 2\ 



qm(m-l)+2j 



(^2„.+2^^2)^(^2„.+2^^2)^ 



n=0 



\rn+j+k=n 



{q^,q^)j{q^,q^)k 
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Then we obtain the exphcit representation 



, m+j+fe=n 



m{m-l)+2jv h^^ ' 



k m 
-X 



m=0 



q—n{l+iy) \ ^m(m— 1) 



^ ^m(m-l) ^ 



{q'"^+',q').{q'^+\q')^_ 



m=0 



j=0 



{q^.q^)Aq^.q^)n- 



X . 



■m-j 



the following identities 



(g2(m-n)^^2) 



Zr?_«i0-l)-2(n-m)j 



leads to 

which gives 

Pnix) = q~ 



2m+2 2\ 

( 2m+2 2n -.yj V^ ^ J n-m jij-l)-2(n-m)j-2mj 

{q ,q)n-m-j-{ i-) 2n q2^ ^ 



( 2m+2 2\ ( 2m+2 2\ ((j2{m-n) 2\ 

yq ^q )n-m-j _ Ig )n-m ' ^ „-2m.? 



■m-j 



n („2m+2 2\ 
■n(l+i/) \ ^ m{m-l)}2. J n-m 

m=0 ''"-"1 



'n-m ('„2?Tt+2 ^2^ ^^2(m-ra) ^2' 



7=0 



j_^2{u-m)j 1 



-n(l+i^) ^ ^ ^m(m— 1) 
m=0 



= g-n(l+i/) ^ ^m(m-l) 



m=0 



n 
m 



n 
m 



n-m ^^2m+2^^2^ ^^2(m-n)^^2^ 



i„2{u-m)j I 



„2m+2 „2(m— n) 



2?^1 



-2n 



g2.^2(.-m) 



Theorem 5 T/te Stieltjes transform of the orthogonality measure /i for the orthogonality poly- 
nomials Pn is given for all z ^ Supp(/x) by 

f djijt) _ 1 
kz-t q(''+^)ju{qV^,q^) 

juiVz, q 



2^ (9-''^-l)(^-\aoo(g'^,g')oo 



g^^'j-.lgl- V^, g2)^(g2^^ ^2)^ _ a^j^^q^^ g2)(g2.^-l^ ^2)^(^2-2.^^ ^2) 

cr,, = 
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Proof. The Stieltjes transform can be obtained from [1, thm. 2.4] 

/ 



Z — t n-^oo Pn{z) 

Using the identity a) in Corollary 3 
g"(^+i)P„(-A2) 

and the identity e) 

g"('^+'-'P„,(-A2) _ 1 
j^(g-"A,g2) - q-2u _ I 

we obtain 



i^(g-"A,g2) 



1 



Qn(-A^) ^ 



Ji/(A,g ) - 



g"(^+i)P„(-A2) 



>(A,g')- 



g(-+i)>(gA,g2) 
Hence, when we use (3) we obtain 

1- Qn(-A^) _ 1 ... 2x 

P„(-A2) - g(-+l)j,(gA,Q2^ ^ 
Now, if A G M5 and using Proposition 1 we have 



q^-j-.{q^--\q^)-jM\Q^) (f^" '^^jjf^^^ 



-2v 



which implies 



If A G g'^Mo then 



Hence we obtain 



lim 

n— >-oo 



lim 



Qn(-A2) i,(A,g2 



+00, 



n^ooP„(-A2) g(-+l)j,(gA,g2 

o„,i-.(cz-"-^A,g2 



lim 

n— ^00 



lim 



Qn(-A2) 



>(9-"A,g2) 



ji/(A,( 



0. 



,-2i/ 



g2'^i_^(gl-'^A,g2)_ 



n^ooP„(-A2) g(-+l)>(gA,g2) 

To complete the proof (case z = A^) we replace A by iX. m 

Corollary 5 The orthogonality measure jx is supported on a denumerable discrete set 

{tk = q~\fk)\ k e N} 

where are the positive simple zero of the function ju{x, q^) numbered increasing. The positive 
mass of fi at t^ is given by 

- 1) 



Ak = 



q^'^j-viq-'jlq^) 
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Proof. From Theorem 5 we deuce that the orthogonahty measure ^ is supported on a denu- 
merable discrete set of the zeros {t^} of jv{q\/z,q^)- So let /x have mass Ak at the point {tk} 
then we have for ah z ^t}. [11, (5.15)] 

^ A, _ 1 

'..r- 2^ jq-^" -l){z-\q^U{q^z,q^)^ 

>^V^, q ) q2uj_^^ql-u^^ q^){z-\q^Uiq^z, q^U ' cTuJAqV^, q^)iq^-z-\q^Uiq^-^-z, q^)^ 

The zeros of jv{q\/z-, q^) correspond precisely to the poles of the left hand side and is the 
residue at this point. ■ 

Remark 3 We can generalize our method to study other families of orthogonal polynomials. 
For instance let V{t) = —^-^ — in (9), then the corresponding orthogonal polynomials Pn{x) 
satisfy 

anPn+l{x) + hnPn{x) + an-lPn-l{x) = xPn{x), 

where 

The family of orthogonal polynomial 

Rn{x) = q^-^'Pn (g-^i+'^^x) 
satisfying the following recurrence relation 

Rn+l{x) + q^^^-^^Rn-lix) = xRn{x) 

with initial conditions 

R-i{x) = 0, Ro{x) = 1. 

The polynomials Rn{x) have the remarkable properties that Rn{q^^x) are also orthogonal and 

are the only modulo a constant (see the reference [2] for more detail). Obviously, using the same 
idea exposed above we can find the corresponding orthogonality measure and generating function 
and probably some other properties. 
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